Berry phase in a non-isolated system 
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We investigate the effect of the environment on a Berry phase measurement involving a spin-half. 
We model the spin+environment using a biased spin-boson Hamiltonian with a time-dependent 
magnetic field. We find that, contrary to naive expectations, the Berry phase acquired by the spin 
can be observed, but only on timescales which are neither too short nor very long. However this 
Berry phase is not the same as for the isolated spin-half. It does not have a simple geometric 
interpretation in terms of the adiabatic evolution of either bare spin-states or the dressed spin- 
resonances that remain once we have traced out the environment. This result is crucial for proposed 
Berry phase measurements in superconducting nanocircuits as dissipation there is known to be 
significant. 
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It was recently suggested Q that it should be possi- 
ble to observe the Berry phase (BP) i 2j in a supercon- 
ducting nanostructure, and possibly use it to control the 
evolution of the quantum state [3, 0- This intriguing 
sugcstion however did not consider the coupling to the 
environment, which is never negligible in such structures 
To truly understand the feasibility of the proposed 
experiment, we must know the effect of the environment 
on the BP. Originally the BP was defined for systems 
whose states were separated by finite energy gaps. Here 
we ask whether a BP can be observed in a system whose 
spectrum is continuous because it is not completely iso- 
lated from its environment. All real systems are coupled, 
at least weakly, to their environment and as a result never 
have a truly discrete energy level spectrum. The usual 
requirement for adiabaticity is that the parameters of the 
Hamiltonian are varied slowly compared to the gap in the 
spectrum. Here there is no gap so naively one would say 
that adiabaticity is impossible and hence the BP could 
never be observed. However experiments have observed 
the BP, both directly and indirectly , so this argument 
must be too naive. We therefore take a simple model in 
which a quantum system, which when isolated exhibits 
a BP, is coupled to many other quantum degrees of free- 
dom. We then ask two questions. Firstly, under what 
conditions can the BP be observed? Secondly, is the ob- 
served BP the same as that of the isolated system? While 
others have investigated systems with a BP coupled to 
other degrees of freedom 0, H, , we believe we arc the 
first to explicitly address these two questions. 

We distinguish between the system and the environ- 
ment in the following way. We have complete experimen- 
tal control over the system, but almost no control over 
the environment. The most that we can do to the envi- 
ronment is to ensure the "universe" (system + environ- 
ment) is in thermal equilibrium, with a temperature T. 
We will assume we have enough control over T to take it 
to zero, and thus prepare the universe in its ground state. 
However any procedure to measure a BP in an isolated 



system must involve measuring a phase difference from a 
superposition of two states. When the system is not iso- 
lated most such procedures involve the mixing of a large 
number of eigenstates (of the universe) , this leads to the 
effects that we discuss below [loj. 

We choose to investigate a spin-half which is coupled 
to both a magnetic field and an environment (a bath of 
harmonic oscillators). Our model is a biased spin-boson 
model 0| with a time-dependent field. When isolated 
from its environment, the spin exhibits a BP if we slowly 
rotate the magnetic field around a closed loop. This 
model, chosen primarily for its simplicity, is extremely 
relevant to a recent proposal for observing a BP in a 
superconducting nanocircuit 0. While we make no at- 
tempt to accurately model the true coupling between the 
nanocircuit and its environment, we believe our results 
give an excellent indication of what to expect in the real 
system. Our work will also be very relevant to realisa- 
tions of the BP quantum computers proposed in Q. 

In this Letter we concentrate on an Ohmic environment 
with the universe initially at zero-temperature |l2| . 
We find that the spin-environment coupling causes the 
spin-eigenstates to become spin-resonances which have 
the following properties, (i) The energy distance be- 
tween them is Lamb shifted by SE. (ii) The higher en- 
ergy resonance exponentially decays to the lower one on a 
time-scale, Ti, and observables containing phase informa- 
tion exponential decay on a timescale T2. (iii) There are 
adiabatic phase-shifts, which divide into two catagories 
with different symmetries; the phase which vanishes when 
the Hamiltonian is time- independent we call <5$bp ! while 
those phase-shifts (and amplitudes) which do not van- 
ish we schematically refer to as $ s hift- The former scales 
with the winding number of the BP experiment, while 
the latter does not (see below). All of these effects go 
like the second power of the spin-environment coupling, 
see eqs. (Jljh®. 

Effect (ii) means that one cannot perform an arbitrar- 
ily long experiment to measure a phase: so we must find 
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the BP from an experiment where the system's Hamil- 
tonian is taken round a closed loop in a finite time pe- 
riod, t p < T 2 . In such an experiment there is typically 
a non-zero amplitude for returning to the initial state 
and this amplitude has a phase. We interpret the lat- 
ter as the sum of a dynamic phase which scales linearly 
with t p , an adiabatic phase (^bp + ^shift) which is inde- 
pendent of t p , and non-adiabatic contributions which are 
proportional to A-t p to some negative power [l^j. Here A 
is the energy difference between the spin-resonances (we 
set U = 1). Thus the BP is present for arbitrary t p , it is 
simply masked by the non-adiabatic contributions unless 
t p is long enough. For the BP to be observed we must 
choose a value for t p which is neither too short nor very 
long, so that it obeys A" 1 <C t p < T 2 . However we then 
actually observe a combination of <1>bp and 'i'shift- To 
distinguish between these two effects we note that when 
we do not rotate the Hamiltonian <1>bp = while $ s hift 
is unchanged. 

Now we ask if the environment's effect on the BP is 
observable. To do this we must first decide what BP we 
would naively expect to observe. There are two possible 
cases to consider: (i) The system evolves in a magnetic 
field that we directly control, then we would expect the 
BP to be given by the solid-angle enclosed by that field, 
$gp. The deviation from this expectation is given by 
(5<&bp in eq. JHJ- For this deviation to be observable it 
must be much larger than the non-adiabatic corrections 
at t p < T 2 ; this means that A • T 2 • <5$ BP > 1- The 
functional form of T 2 and <5$bp, in 10} and (JSJ, have 
the same dependence on the strength of the coupling to 
the environment, C. Thus the condition reduces to one 
dominated by the dependence on 7, where 7 (defined be- 
low eq. |J5J) characterises the environment. We conclude 
that there is a wide range of values of 7 for which we 
can observe (5<f> B p. (ii) The second case is more compli- 
cated, but is relevant to the superconducting nanocircuit 
in . There we have no independent measure of the bare 
spin Hamiltonian, the control parameters (gate voltages 
and magnetic fluxes) enter the spin Hamiltonian in com- 
bination with unknown constants (capacitances and in- 
ductances). Thus we know nothing about the bare spin- 
eigenstates, or the solid angle that they enclose when 
we vary the experimental parameters. However we can 
measure the spin resonances in the presence of the en- 
vironment as a function of the experimental parameters. 
Then one might predict the observed BP is given by the 
solid-angle enclosed by these spin-resonances. This pre- 
diction is given above JJJ; it is of a similar form to the 
correct result, but contains a very different function of 
the distribution of oscillators in the environment. The 
deviation from this expectation is given by SQ'^p, for it 
to be observable we require that A-T 2 -5& BP S> 1. Again 
this reduces to a function independent of C, where <5$ BP 
is observable over a wide range of 7. Finally, we assume 
we measure $ s hift when for a time-independent Hamilto- 
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FIG. 1: Evolution in step (b) of the experiment, in Lab. and 
rotating frames. The primed-basis and plus-basis are both 
shown. 

nian before carrying out the BP experiment. Then we do 
not require <5$ B p (or <5<i>gp) to be larger than $ s hift for 
it to be observable. 

To be concrete we assume here that the BP is mea- 
sured using the spin-echo method[j|. We consider an 
experiment where we start with the field along the z-axis 
and the universe (spin+oscillators) is in its ground state 
[l4|. (a) The field is then (instantaneously) prepared at 
its initial value, Bq, which is at angle 9 to the z-axis, 
At the same time the spin is (intantaneously) placed in 
the state -^j(IT) + ID) relative to Bo- Then (b) we adia- 
batically rotate the magnetic field, B(t), n times around 
a closed loop with constant angular velocity, uj = zu, 
(see Fig. for a time period, t p = 2nn/ui. We call n 
the winding number. After which (c) the spin is flipped 
and (d) the field is rotated with angular velocity — w for 
time t p . Finally (e) the spin is flipped again and (f) 
the spin state is measured. By "flip the spin" we mean 
IT) IDj wnere the f and J, are relative to the direc- 
tion of the B-field at that time. This can be achieved 
by applying a instantaneous 7r-pulse oriented along the 
y-axis. By instantaneous we mean much faster than the 
fastest oscillator in the environment. We ask what the 
probability is that the final spin-state, after carrying out 
(a)-(f), is in a given direction in the plane perpendicular 
to Bq. For an isolated spin-half, the probability of the 
final spin-state being ^(e i4/2 |T) + e" i?/2 | j)) is 



P(0 = \ 



cos(£-4< p 



(1) 



where all spin-states are defined relative to the axis of 
the field Bq. Measuring this probability as a function 
of £ yields the BP for an isolated spin, <f> 
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cos 9) . We wish to know what we observe if we carry 
out the same measurement for a spin which has been 
weakly coupled to a bath of oscillators throughout the 
experiment. 

The Hamiltonian we consider contains the spin-half in 
the above time-dependent magnetic field, B(t), which 
is also coupled to a bath of harmonic oscillators with 
frequencies {fij}. Writing it in terms of creation, lv, and 



3 



annihilation, 6, operators for the oscillators, 

H(t) = -|B( t ).a + ^%(St a S, > + i 



9 \ ^ Cq 

2 (2m%) 



T/2 (i + *° , ( 2 ) 



where j is summed over all oscillators and a is summed 
over the (x, y, z) components of the oscillator. The 
number of oscillators with frequency fl to Q + dfi is 
p(il)dil. The spectral density [ll| is given by J(f2) = 
Ej 7i"(5C) 2 (2m%) _1 <*(fi - Oj) = 7r( 5 C) 2 p(fi) (2mO) _1 . 
Here we restrict ourselves to z-axis spin-environment cou- 
pling with C a = CS az 0. Then for B ■ a = £<j z , 
the exact ground state of the universe 14] is simply 
T}llj|0 J T ) where oscillator j is in the ground state, 
|0i), of the harmonic potential centred at (0, 0,^gC). 
We consider an Ohmic bath of oscillators with J(f2) = 
T;C 2 ttex-p [— f2/f2 m ], and work in the limit of small di- 
mensionless coupling C = gC{a/ra) 1 / 2 <C 1. 

The time-dependence in @ makes the problem un- 
pleasant, however we remove this by going to the primed- 
basis which rotates with the £?-field. In this non-inertial 
basis the spin experiences an effective field [Bq + g^ 1 ^) ■ 
For our problem the effective field is B + for < t < t p 
(shown in Fig. QJ, and B^ for t p < t < 2t p , where B± — 
(Bo ± g~ x ijjz). Having removed the time-dependent, we 
calculate the evolution of the system in the frame which 
has its z-axis parallel to the field (either B + or B_) these 
frames we call the plus- and minus-basis respectively (the 
former is shown in Fig. ^) . Finally we rotate back to the 
lab-frame to evaluate observables. 

Before we give a detailled explanation of how we cal- 
culate the spin's evolution in the presence of the envi- 
ronment, we giving our results. The anisotropic nature 
of the coupling results in P(£) containing 0[(7 2 ]-terms 
which go like exp[±i<?_Bi p ]. To simplify the resulting ex- 
pressions we average t p over a range > (gB)^ 1 to remove 
these terms, then 
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where $bp = 3> 



|ks| e" 2tp/T2 cos (f + 4$ B p - arg[ K2 ]) 

\k 3 \ (2e- 2 *p/ T2 - e ~ 4tp/T2 ) cos (£ - arg [n 3 ]) 

(3) 

<5$bp- For compactness we have 
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dropped an uninteresting real 0[C 2 t p ) term from the first 
exponent while retaining such terms elsewhere. The ks 
(which were schematically refered to as $ s hift above) are 
0[C 2 } and so are comparable to <5$bp; however they are 
independent of the time-dependence of H(t), and hence 
independent of the winding number, n. We find, 

Ts" 1 = (2T 1 )- 1 = fC ,2 O m7 e- 7 sin 2 6l (4) 
6*bp = \nC 2 [/'( 7 )-2 7 - 1 /(7)]sin 2 0cos0 , (5) 



where 7 = gB/£l m . The function f(x) = xe x Ei(— x) + 
xe _x Ei(x) where we define Ei(x) as the principal-value 
of the Exponential integral, dte^/t, and f'(x) = 
df(x)/dx. Eq. (J5J is simply the w-dependent term in 
the Lamb shift of the energy when in the rotating frame. 
This generates a term of 0[n-i p ] in the phase which in 
the laboratory frame is a contribution to the BP. 
The n-independent factors are 

«i = ±C ,2 7re- 7 sin0 [cos0 + ±sin0] 

K2 = ^C 2 ( 7 -V(7)+i^)sin 2 

« 3 = i(7 2 [( 7 - 1 /( 7 )+i7re- 7 )sin0cos0-2 7 - 1 sin0] 



K4 



\C 2 [e 7 Ei(- 7 )sin0cos( 
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Now we check that the BP is not simply given by the 
solid-angle enclosed by the spin-resonances, //this were 
the case then the BP for this experiment would be 
fnC' 2 7 _1 /( 7 ) sin 2 0cos0, the correct result de- 
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viates from this prediction by 



<S$ BP = fnC' 2 /'( 7 )sin 2 0cos0 , 



(7) 



for most 7 this deviation is significant. 

We now discuss the method we use to obtain these 
results. The Hamiltonian in the primed-basis is time- 
independent and is given by 



n'(B ± ) 
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If we write the spin's initial density matrix as p and 
the oscillators initial density matrix as pg sc , then we 
are interested in the spin density matrix at time t, 
after we have traced over the oscillator states, p t = 
tr OS c^t [po ® Po sc ) i where U t is the evolution opera- 
tor. We find it helpful to write the spin density ma- 
trix as a vector p whose elements are (pn, P12, P21, P2i)- 
Then the spin evolution equation (after the oscillators 
have been traced over) can be written as p t = K(£)p , 
where this defines K(f) as a four-by-four matrix which 
gives the time evolution of the elements of the spin's 
density matrix. The initial state of the oscillators en- 
ters in the functional form of the elements of K(f). For 
the experiment described above eq. JIJ we need to cal- 
culate p 2tp = K fli P K(B + ,t p ) K fli P K(B_,t p ) p . The 
spin-flip is assumed to be fast enough to leave all the os- 
cillators unchanged while flipping the spin, then in the 
primed-basis K fllp simply has "l"s on the off-diagonal 
and "0"s elsewhere. This leaves the calculation of the 
propagation matrix K(_B + ,< p ), we can find K(-B_,t p ) 
by reversing the sign of uj throughout. For weak cou- 
pling to the bath it is natural to work in the plus-basis 
(see Fig. UJ, which has its z-axis parallel to B + , in this 
basis ~K(B + . t p ) becomes diagonal if C — ► 0. Finally the 
coupling between spin and oscillators in the plus-basis 
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FIG. 2: The four classes of O [C 2 ] -contributions to E+(r) 
are shown here. The upper (lower) line is the Retarded (Ad- 
vanced) non-interacting (C a = 0) spin-propagator. Every- 
thing is in the basis where these propagators are diagonal (the 
plus-basis for E + ). The spin-state is fj, with j (J.) = +1(— 1). 
The dotted-lines are the spin-spin interactions after we have 
traced out the oscillators which couple to the spin via a a . 

at time t is C+(t) = C'(i)R + , where R+ is the SO(3) 
rotation from the primed-basis to the plus-basis. 

Now we use the real-time transport method |16| to 
write the following differential equation for K + (t), 



d t K+(t) = -iE+K+(i) + f drS+(T)K+(/j 

Jo 



r) , (9) 



where all bold symbols are 4x4 matrices. The ma- 
trix E + gives the evolution of the propagation matrix 
when there is no coupling to the bath. Because we are 
in the plus-basis it is diagonal with E^ — E^ = 
and E+ 2 = -E+ 3 = -gB+ = -g\B + g~ l u\. The 
matrix S + (t) is the contribution of all irreducible di- 
agrams with one or more interactions with the bath of 
oscillators. Equation © is exact, however to proceed 
we treat this equation to first order in C 2 . Thus in the 
integral on the left hand side of © we treat S + (t) to 
first order in C 2 and K + (t — r) to zeroth order. So we 
can write K + (i — r) ~ K + (i)Kj(— t) where the cor- 
rections to the approximations are 0[C 2 ] and so can be 
ignored. Now X + , which is evaluated below, is dom- 
inated by small-r so we take the upper limit on the 
integral to infinity. The error we make in doing so is 
[(Onip)" 1 ] which we neglect. This systematic ap- 
proximation results in the interaction becoming local in 
time. Then we get d t K+(t) = [-iE+ + X+] K+(t) where 
X + = J °° drS + (T)K,|(— r), and diagonalise the matrix 
(-iE+ + X+) to find K+(t). 

Now we briefly discuss the evaluation of S + to lowest 
order in C 2 . At this order we need only consider irre- 
ducible diagrams with a single interaction with the os- 
cillators. When the oscillators are traced out they leave 
an interaction between the spin at time t and time t — r. 
The resulting first-order irreducible diagrams are shown 
in Fig. |21 The contribution to S + of the diagram with 
an interaction via a a at time t and another via a a i at 
time t' after we have summed over the Ohmic bath, is 



weak enough that gB ^> T 2 
4°p + <^bp 




m \ (10) 



where re = ±1, % = ±1, and other variables are shown in 



Fig. The upper (lower) term in {• ■ •} is applicable if 
the relevant vertex is on R (A), k is +1 (—1) when the 
a' vertex is on R (A), x 1S +1 (~ 1) if the interaction is 
R-R or A-A (R-A or A-R). 

In conclusion, the BP can be observed in a non- 
isolated system, if the coupling to the environment is 

The adiabatic phase is 

$ s hift, but <& s hift is not considered a BP 
because it does not vanish when n — 0. So the BP dif- 
fers from that of a isolated spin by 5$bp, given in Eq. 
(0). The proportionality of ^^bp to n hints that it has 
some geometric character, however it is a function of the 
environment's spectrum and thus the total BP is not a 
simple geometric quantity. 
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